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Abstract—A statistical analysis of the single-delay multiplica- 
tion based frequency estimator is treated here. The single-delay 
multiplication frequency estimator for complex single sinusoid 
signals is a phase averaged estimator, which is similar to the Kay’s 
estimator. Here we provide a study on the statistical distribution 
of the frequency estimates made by the estimator, and verify the 
analytical results using simulations. It is shown that the analytical 
results hold true even at very low signal to noise ratio levels. 
In deriving the probability density function of the frequency 
estimates, we explore the cross product noise terms resulting 
due to the multiplication operation of the estimator by making 
valid assumptions. 


I. INTRODUCTION 


The estimation of frequency of a single tone complex signal 
is a well treated problem in literature [1-5]. Methods vary 
from open loop techniques to robust closed loop or feedback 
techniques. The most attractive estimator in terms of jitter 
performance is the maximum likelihood (ML) based frequency 
estimator, which involves high complexity in terms of practical 
implementation [3-5]. The single-delay multiplication (SD- 
M) based frequency estimator on the other hand is less 
efficient with respect to the jitter performance of the ML 
frequency estimator but with reduced complexity. The single- 
delay multiplication estimator is a kind of phase averaging 
estimator, which is similar to the Kay’s estimator [1]. Kay, 
in his paper, provided this technique also with some other 
modified versions of it by incorporating weighting functions. 
In this paper we analyse the proposed estimator by performing 
statistical analysis on it. The statistical distribution of the 
frequency estimates at the output of the estimator is derived. 


Il. SINGLE-DELAY MULTIPLICATION BASED 
FREQUENCY ESTIMATOR 


The formation of the estimator is straight forward. Consider 
the received complex sinusoidal signal corrupted with addic- 
tive complex white Gaussian noise, 


r[n] = Aexp{jQn + jo} +7 (1) 
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Fig. 1. 
estimator. 


Block diagram of the single-delay multiplication based frequency 


where, A, 2 and @ are all deterministic constants, and 77 is the 
complex white Gaussian noise process given by, 


N= Net jns (2) 


where, n. and n, are two independent zero mean white 
Gaussian random processes with o? variance and a double 
sided power spectral density of No each. Then the single-delay 
multiplication based frequency estimator is given by, 


L 

Q = arctan {> r[njr*[n — u} (3) 
n=1 

where, L+1 is the number of samples required per estimate. 

For the noiseless case, the estimator in (3) estimates the precise 

discrete frequency in radians given by, 


Q = 2n fT, (4) 


where, 7’, is the sampling period of the discrete signal and 
f is the frequency in Hertz. The arctan function is the four- 
quadrant inverse tan operation whose statistical properties are 
treated later in the paper. Fig-1 shows the block diagram 
of the frequency estimator. Kay’s analysis on the frequency 
estimator, similar to the one given in (3), was based on a 
linearised expression on the additive noise component for high 
signal to noise ratio. The model that was used by Kay in [1] 
is given by, 


r[n] = Aexp{jQn + jO +v} (5) 


where, v is a real zero mean white Gaussian random process. 
This assumption disregards the nonlinear effects caused by the 
arctan operation in the estimator. The effect is not apparent 
at high signal to noise ratio levels, however at low signal to 
noise ration levels the nonlinear characteristics of the estimator 
has to be considered in order to provide an appropriate model. 
Here, in our analysis, we consider the nonlinear characteristics 
of the arctan operation and develop a statistical model which 
holds true even at low signal to noise ratio levels. The 
analytical results are verified using simulations. 


III. NOISE ANALYSIS 


Consider the block diagram of the estimator given in Fig-1. 
Let the signal at the output of the multiplier in the figure be 
z[(n]. Then, z[n] is given by. 


z[n] = r[n]r*[n — 1] = A?[exp{jQ}+Pi+ P+ © 


where, the terms, P,, Pz and € are given by, 


1 
Pi = 5 exp{jQn + jO}n*[n — 1] (7) 
1 
Py = i exp{—jQ(n — 1) — jO}n[n] (8) 
and 1 
€= Fanln|n*[n - 1] cc) 


The processes P, and P2, are two linear transforms of the 
original noise process 77, and therefore they follow a Gaussian 
distribution with the inphase and the quadrature components 
having zero mean and a variance of o7/A?. Let us now 
consider the noise cross product term € given by (9). The 
expression in (9) is expanded and rewritten in terms of its 
inphase and quadrature components as, 


€= lla +] + JI - al} (10) 
where, 
= neln}ngin —1 (11) 
&=n,lnlnele— 1 (12) 
= aint (13) 
4 =n,[n]n,[n —1 (14) 


where, n, and n, are the inphase and the quadrature compo- 
nents of the original noise defined in (2). Here, we need to 
analyse the statistical properties of the terms in (11) and (12). 
Consider the term &), since n,[n] and n,[n — 1] are zero mean 
Gaussian random variables with variance o?, the joint density 
function of n.[n] and n,[n — 1] is given by, 


1 
27r| det K|!/2 


1 
f(A) = exp{ 5 (A7K7!X)} (15) 
where, 


A= [n,[n]_ n-[n — 1]]" (16) 


and, K is the covariance matrix given by, 


K= i ei 


O11 oO 


(17) 


In (17), o11 is the covariance of n,[n] and n,[n — 1], and we 
recognize that, 


on = E[n.[n|n-[n — 1]] = R,(Ts) (18) 


where, R,,(7) is the autocorrelation function of n, at T = Ty. 
Assuming that the noise process is ideally low passed filtered 
with a rectangular brick wall type of filter with a bandwidth 
of B; Hz, then the autocorrelation function R,(7) can be 
expressed as, 


R(T) = 2NoBisinc(2B;r) (19) 


The autocorrelation function in (19) becomes zero when the 
filter bandwidth B; satisfies the condition, 
2=— 

2Ts 
where, m is a positive integer. For discrete systems, when 
signals are ideally sampled at a sampling frequency of 1/T; 
Hz, equation (20) holds true with m = 1. Thus, when equation 
(20) is satisfied, the covariance of n,[n] and n.[n—1] becomes 
zero, resulting in a joint distribution function given by, 


2 2 


(20) 


AQ) = sa us { 


Now, by using the result in (21) and making the necessary 
linear transformation on n.[n], or otherwise, we can find the 
probability density function of &1, given by, 


x 1 
fe, (&1) =i 2n| Blo? exp 


where, 3 = n,[n]. A closed form expression for the 
probability density function in (22) is not explicit. Hence, 
numerical techniques were used to solve the integral in (22). 
The results were verified using simulations and are shown in 
Fig-2. Even though a closed form expression for the pdf was 
not obtained in (22), the mean and the variance of £, are easily 
obtained by direct analysis. The variable €; has zero mean 
because of the symmetrical nature of the pdf in (22). The 
variance is computed by solving a double integral involving 
(22), and is given by o*/A*. With these first order statistics of 
€,, we can determine the same for the real and the imaginary 
components of € summed over L terms, given in (10). Let, 


x= ol + & 


n=1 
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\ dB (22) 


(23) 


Note that €; and & are intrinsically functions of n as 
described in (11) and (12) respectively. Since the variances 
of €, and & are bounded, using the central limit theorem [7], 
we find that the distribution of is a Gaussian for large values 
of L (typically for ZL > 5 from simulations). The mean and 
the variance of y are given by m, = 0 and o% = 2Lo*/A*, 


noise variance = 2.5 + simulation 
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respectively for L > 5. The above results are verified by 
simulation and are shown in Fig-3 for L = 10 and L = 60. 
Now using the noise analysis performed in this section, we 
determine the distribution of the frequency estimates made by 
the SD-M estimator in the following section. 


IV. STATISTICAL DISTRIBUTION OF THE 
FREQUENCY ESTIMATES 


The statistical properties of the frequency estimates made 
at the output of the SD-M estimator are of great interest to 
us. Especially one needs to know how the estimates behave 
Statistically. Let us rewrite the expression in (3) taking into 
account the additive noise component. 


L 
Zz =) aln| = A?{Lexp{j} +x +x2+x3} 24) 
n=1 
where, 
L L L 
x1 = SOP .x2=S > Prandxs=S°E (25) 
n=1 n=1 n=1 


In the above equation, x; and x2 are Gaussian distributed ran- 
dom processes, each with means zero and variances Lo? /A?. 
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Fig. 4. PDF of the frequency estimates Q, SD-M estimator 


The real and the imaginary components of 3 have similar 
statistical properties as y. Now, we could rewrite (24) as, 


Z,= 44+ 922 (26) 


In (26), Z; and Z_ are two Gaussian random processes with 
means j1; and ji2 respectively, and variances o3 each, given 
by, 


pi = A? L cos(Q) , wa = A? Lsin(Q) (27) 
and 
2 2 a 


Then, the frequency estimate at the output of the SD-M 
estimator is given by the four-quadrant arctan function of 
vam 


Q= arent?) (29) 
Z\ 

The equation in (29) is treated as a transformation of multiple 

random variables in order to determine the output distribution. 

Since the distributions of Z, and Z 2 are known, we can 

determine the statistical distribution of the frequency estimates 

made at the output of the estimator, which is given by, 


fa(Q) = ferea(21s22)IJ(21,22)| der 30) 


where, 22 = 21 tan((), and |J(z1, 22)| is the determinant of 
the Jacobian matrix for the transformation given in (29), which 
is given by, 


Oz, Oz | 
| J(21, 22)| = . aT (31) 


The probability density function of Q is then given by solving 
the integral in (30) and the results are given in (32), (33) 
and (34). Since the trigonometric arctan function has got two 
discontinuities within the range of —7 to 7 , one at —7/2 and 
the other at 7/2, the integral in (30) needs to be evaluated for 
three different regions. The analytical results obtained for the 
statistical distribution of the frequency estimate at the output 
of the SD-M estimator is verified using simulations, and the 
results are shown in Fig-4. The proposed statistical model for 
the estimator matches very closely to the simulated results. 


exp(-s/2) E exp (4) (4) evm)| 


fa(O) = { exp(—s/2)|  +exp (4) (A) (Q-@(VA))] tor - 2/2 << 4/2 


exp(—</2) E exp (48) (42) @ (VD) | 


is 2 
stad «lan Rn 


for t/2<QA<n 


for —17 <Q <—1/2 


[tan (« + 0) ba + p12] : 


2A= A, = 


where ,s = . - 
2 [1 + tan? (2)| o3 


[tan (-7 + 0) bai + p12] : 


Ag = and ,Q(a) = 2 ie exp (—u?/2) du 


[1 + tan? (-7 + 0)| o3 


V. CONCLUSION 


Statistical analysis of the single-delay multiplication based 
frequency estimator was treated here. The distribution of 
the frequency estimates made at the output of the estimator 
was determined and verified using simulations. The statistical 
model is very precise even at very low signal to noise ratio 
levels, and is of great use in receiver designs to estimate the 
jitter and the behavior of the estimator when operating under 
low signal to noise ratio conditions. 


ACKNOWLEDGMENT 


The authors would like to thank the Wireless Signal Process- 
ing program at the National ICT Australia. National ICT Aus- 
tralia is funded through the Australian Government’s Backing 
Australia’s Ability initiative, in part through the Australian 
Research Council. Kandeepan is also an adjunct fellow with 
the Research School of Information Sciences and Engineering 
at the Australian National University, and would like to thank 
the school for the extended support. The authors would also 
like to thank Dr. Leif Hanlen from the National ICT Australia 
for his assistance in preparing the paper. 


REFERENCES 


[1] S. Kay, “A Fast and Accurate Single Frequency Estimator,’ JEEE Trans- 

actions on Acoustics, Speech and Signal Proc, vol. 37 No.12, pp. 1987- 

1990, Dec 1989 

[2] S. Kay, “Fundamentals of statistical signal processing,estimation theory.” 

Englewood Cliffs, N.J: Prentice-Hall, 1993. 

[3] D. C. Rife and R. R. Boostyn, “Single-tone parameter estimation from 

discrete-time observations”, JEEE Trans on Inform Theory, Vol IT,No 

5,pp 591-598, Sep 1974. 

[4] U. Mengali and A. N. D’Andrea, Synchronization Techniques for Digital 

Receivers, Plenum Press, 1997. 

[5] H. Meyr, M. Moeneclaey and S. A. Fechtel, Digital Communication 

receivers, Synchronisation, Channel Estimation and Signal Processing, 

John Wiley and Sons, 1998. 

[6] S. Kandeepan, “Synchronisation Techniques for Digital Modems”, Phd 

thesis, Uni of Tech, Sydney, July 2003 

[7] Anthanasios Papoulis, Probability, Random Variables, and Stochastic 
Processes, McGraw-Hill, 2nd Ed, New York 1984. 


[1 + tan? (= + 0)| ge 


(32) 


(33) 


(34) 


